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ABSTRACT 

We present the luminosity function of 90/^m selected galaxies from the Eu- 
ropean Large Area ISO Survey (ELAIS), extending to z = 0.3. Their lumi- 
nosities are in the range 10 9 < h^L/L e < 10 12 , i.e. non-ultraluminous. 
From our sample of 37 reliably detected galaxies in the ELAIS SI region 
from the Efstathiou et al. (2000) Sqo > lOOmJy database, we found opti- 
cal, 15/im or 1.4GHz identifications for 24 (65%). We have obtained 2dF 
and UK Schmidt FLAIR spectroscopy of 89% of IDs to rigid multivariate 
flux limits. We construct a luminosity function assuming (a) our spectro- 
scopic subset is an unbiased sparse sample, and (b) there are no galaxies 
which would not be represented in our spectroscopic sample at any red- 
shift. We argue that we can be confident of both assumptions. We find 
the luminosity function is well-described by the local 100/im luminosity 
function of Rowan-Robinson, Helou & Walker (1987). Assuming this lo- 
cal normalisation, we derive luminosity evolution of (1 + £) 2 - 45 ±°- 85 (95% 
confidence). We argue that star formation dominates the bolometric lu- 
minosities of these galaxies and we derive comoving star formation rates 
in broad agreement with the Flores et al. (1999) and Rowan-Robinson et 
al. (1997) mid-IR-based estimates. 

Key words: cosmology: observations - galaxies: formation - infrared: 
galaxies - surveys - galaxies: evolution - galaxies: star-burst 



1 INTRODUCTION 

The study of the star formation history of the Uni- 
verse is an extremely active field, in which much of 
the current debate is centred on the uncertainties in 
dust obscuration. Selecting star forming galaxies in 
the UV is relatively cheap in observing time, but 
such selection makes the samples extremely sensi- 
tive to dust obscuration. Detecting the reprocessed 



starlight requires sub-mm surveys (e.g. Smail, Ivison 
& Blain 1997; Hughes et al. 1998; Barger et al. 1998, 
1999; Eales et al. 1999, Peacock et al. 2000; Ivison et 
al. 2000a,b, Dunne et al. 2000) and/or space-based 
mid-far-IR surveys (e.g. Rowan- Robinson et al. 1997, 
Taniguchi et al. 1997, Kawara et al. 1998, Flores et 
al. 1999, Puget et al. 1999, Oliver et al. 2000; see 
e.g. Oliver 2000 for a review). As such, the European 
Large Area ISO Survey (ELAIS) is well-placed for 



© 2000 RAS 



2 Stephen Serjeant, et al 



the study of the evolution and obscuration of the star 
formation in the Universe, and strong constraints at 
z S 1 are possible from ELAIS (Oliver et al. 2000). 

The scientific aims and strategy of ELAIS are 
presented in detail in Paper I (Oliver et al. 2000). In 
Papers II and III (Serjeant et al. 2000, Efstathiou et 
al. 2000) we presented respectively the ELAIS Prelim- 
inary Analysis source counts from the CAM (6.7/im 
and 15/xm) and PHOT (90^im) instruments on ISO. 
The 90/xm sample covered 11.6 square degrees, and 
the source counts were found to agree well at the 
bright end with the IRAS 100/im counts. Excellent 
agreement was also achieved with a parallel indepen- 
dent pipeline (Surace et al. in prep.). 

In this paper we present the first 90/mi luminos- 
ity function from our initial spectroscopic campaigns 
in the ELAIS SI field (Oliver et al. 2000 in prep., 
Gruppioni et al. 2000 in prep., La Franca et al. 2000 
in prep., Linden- V0rnle et al. 2000 in prep.). Section ^ 
defines our sample. In section ^| we use this catalogue 
to derive a 90^tm luminosity function. We discuss the 
implications of our results in section ^| We assume 
H = 65 km s _1 Mpc -1 , Q = 1, A = throughout. 



2 SAMPLE SELECTION AND DATA 
ACQUISITION 

The parent sample for this study is the Preliminary 
Analysis catalogue of Efstathiou et al. 2000 (Paper 
III), with 90/im fluxes satisfying Sgo > 100m Jy. The 
completeness of this sample falls approaching this 
limit but has been well-quantified with simulations 
(figure 6 of Paper III). We will use the calibration 
adopted in Paper III. The 90/mi flux calibration is 
still uncertain to within ~ 30%, as discussed in Pa- 
per III; however, in order to be consistent with the 
100/im IRAS source counts (assuming Sgo ^ Sioo) 
the flux calibration has to be very close to the Paper 
III value. Therefore when we compare our 90/^m lu- 
minosity function to the local 100/im luminosity func- 
tion below, we will treat the flux calibration as being 
accurately known. Nevertheless, in predicting the red- 
shift distributions of other 90/im surveys, we will not 
neglect the flux calibration systematics. 

We restrict our study to the ELAIS SI field with 
90/im, 1.4GHz and 15/im survey coverage, which cov- 
ers 3.96 square degrees and has the most extensive 
available optical spectroscopy. Large-scale-structure 
variations are negligible in the comoving volume of 
the total ELAIS survey, but the ELAIS SI field has 
slightly lower source counts than the global average. 
Consequently we reduced the effective area by 20% 
(a factor of 0.8) to make the SI counts normalisation 
match the whole survey, to account for any possible 
large scale structure fluctuations. (This correction is 
smaller than the errors in the luminosity function de- 
rived below.) We sought identifications of this sample 
using (a) the APM (Automatic Plate Measuring ma- 
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Figure 1. Luminosity - rcdshift plane. The 100m Jy flux 
limit is marked. 



chine) , (b) sub-mjy decimetric radio sources from our 
Australia Telescope Compact Array (ATCA) imag- 
ing (Gruppioni et al. 1999), (c) our ELAIS ISOCAM 
6.7/^m and 15^tm catalogues (Serjeant et al. 2000, Pa- 
per II). By randomising our 37 catalogue positions 
in the cross-correlation, we found that < 1 optical 
APM identifications brighter than R = 17 should 
occur by chance within the error ellipses of all our 
PHOT sources, and similarly restricted the associa- 
tion radius for the CAM and radio catalogues. Optical 
spectroscopy was obtained at the 2 degree field (2dF) 
at the Anglo- Australian Telescope (AAT), the Fibre 
Linked Array Image Reformatter (FLAIR) spectro- 
graph at the UK Schmidt and with the Danish Faint 
Object Spectrograph and Camera (DFOSC) at the 
Danish 1.54m at ESO. The optical spectroscopic cat- 
alogue will appear in future papers in this series. 
Of the 37 reliably detected galaxies, 24 had multi- 
wavelength identifications of which 18 are above the 
combined multivariate flux limits discussed above. Of 
these we have obtained optical spectroscopy of 89%. 
The ELAIS SI 90/im sample is presented in table [l] 
with the cross-identifications and available redshifts. 
This table also lists redshifts for four identifications 
for which the identification is uncertain. These galax- 
ies fail the strict selection function defined below and 
are excluded from the luminosity function analysis 
below. Two galaxies (ELAISP90_J00372 1-434228 and 
ELAISP90_J003021-423657) are Seyfert II galaxies 
and two further galaxies (ELAISP90_J003431-433806 
and ELAISP90_J003133-425100) have an early type 
spectrum; all the remainder (including the four ex- 
cluded identifications) show starburst emission line 
spectra. 



3 RESULTS: THE 90 MICRON 
LUMINOSITY FUNCTION 

Our selection function is as follows. In order to be in 
the spectroscopic target list, a 90/im galaxy must ei- 
ther have (a) a reliable optical identification brighter 
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Name RA Dec £go/Jy "Sis/mjy Si^/mJy R z z max L/Lq 

(J2000) (J2000) 



ELAISP90.J 003857-424358 


00 


38 


57.7 


-42 


43 


58.8 


0.43 


2.9 


0.964 














ELAISP90.J003056-441633 


00 


30 


57.0 


-44 


16 


33.6 


0.35 


4.9 




13.7 


0.019 


0.0351 p 


2.8 


X 


10 9 


ELAISP90.J003510-435906 


00 


35 


10.5 


-43 


59 


6.0 


0.41 


9.6 




12.1 


0.136 


0.272 p 


1.7 


X 


10 11 


ELAISP90.J003431-432613 


00 


34 


31.6 


-13 


26 


13.2 


0.59 






14.0 


0.053 


0.127 p 


3.7 


X 


10 10 


ELAISP90.J003914-430437 


00 


39 


14.8 


-43 


04 


37.2 


1.79 


16.2 




11.6 


0.012 


0.049 p 


5.7 


X 


10 9 


ELAISP90.J003021-423657 


00 


30 


21.1 


-42 


36 


57.6 


1.17 


22.3 




16.0 


0.149 


0.493 p 


6.0 


X 


10 11 


ELAISP90.J003459-425718 


00 


34 


59.0 


-42 


57 


18.0 


0.55 


5.5 


1.371 


14.9 


0.055 


0.101 p 


3.7 


X 


10 10 


ELAISP90.J003134-424420 


00 


31 


34.5 


-42 


41 


20.4 


0.31 






13.1 


0.027 


0.048 p 


5.1 


X 


10 9 


ELAISP90.J003242-423314 


00 


32 


42.9 


-42 


33 


14.4 


0.24 


2.6 


0.463 


15.2 


0.053 


0.081 p 


1.5 


X 


10 10 


ELAISP90_J003615-424344 


00 


36 


15.4 


-42 


43 


44.4 


0.31 






17.2 


0.115 










ELAISP90.J003516-440448 


00 


35 


16.1 


-44 


04 


48.0 


0.29 


















ELAISP90.J003741-440227 


00 


37 


41.6 


-44 


02 


27.6 


0.21 






19.9 


0.348 










ELAISP90.J003431-433806 


00 


34 


31.2 


-43 


38 


6.0 


0.75 




0.595 


18.4 


0.199 


0.357° 


6.8 


X 


10 11 


ELAISP90.J003405-423816 


00 


34 


5.3 


-42 


38 


16.8 


0.23 


















ELAISP90.J003358-441 102 


00 


33 


58.2 


-44 


11 


2.4 


0.23 


5.0 




18.9 


0.305 


0.446° 


5.2 


X 


10 11 


ELAISP90.J003501-423852 


00 


35 


1.9 


-42 


38 


52.8 


0.28 




0.606 


15.9 


0.054 


0.090 p 


1.8 


X 


10 10 


ELAISP90.J003254-424608 


00 


32 


54.2 


-42 


46 


8.4 


0.46 


2.5 


0.469 


17.7 


0.190 


0.207° 


3.8 


X 


10 11 


ELAISP90.J003531-423314 


00 


35 


31.4 


-42 


33 


14.4 


0.16 


















ELAISP90.J003635-430148 


00 


36 


35.7 


-43 


01 


48.0 


0.12 


4.5 


0.358 


21.5 












ELAISP90_J003019-432006 


00 


30 


19.8 


-43 


20 


6.0 


0.20 


















ELAISP90.J003839-422324 


00 


38 


40.0 


-42 


23 


24.0 


0.14 


















ELAISP90.J003014-440509 


00 


30 


14.3 


-41 


05 


9.6 


0.10 


3.6 
















ELAISP90.J003912-431203 


00 


39 


12.2 


-43 


12 


3.6 


0.15 


5.0 




14.8 












EL AISP90.J003 133-425 100 


00 


31 


33.9 


-42 


51 


0.0 


0.11 


6.3 




18.0 


0.211 


0.222 p 


1.1 


X 


10 11 


ELAISP90.J003304-425212 


00 


33 


4.3 


-42 


52 


12.0 


0.19 


1.4 




15.6 












ELAISP90.J003046-432204 


00 


30 


46.4 


-43 


22 


4.8 


0.11 


6.2 




15.6 


0.073 


0.075 p 


1.3 


X 


10 10 


ELAISP90.J003027-433050 


00 


30 


27.7 


-43 


30 


50.4 


0.11 


4.9 




15.7 


0.072 


0.074 p 


1.2 


X 


10 10 


ELAISP90.J002845-424420 


00 


28 


45.5 


-42 


44 


20.4 


0.11 


















ELAISP90.J003719-421955 


00 


37 


19.4 


-42 


19 


55.2 


0.16 


















ELAISP90.J003725-424554 


00 


37 


25.2 


-42 


45 


54.0 


0.12 


















ELAISP90.J003731-440758 


00 


37 


31.3 


-44 


07 


58.8 


0.32 






20.0 












ELAISP90.J003625-441 127 


00 


36 


25.0 


-44 


11 


27.6 


0.14 




4.265 














ELAISP90.J003721-434228 


00 


37 


21.5 


-43 


42 


28.8 


0.12 


17.3 


0.451 


17.4 


0.225 


0.248 p 


1.1 


X 


10 11 


ELAISP90.J003348-433032 


00 


33 


48.3 


-43 


30 


32.4 


0.15 


















EL AISP90.J003 149-423628 


00 


31 


49.2 


-42 


36 


28.8 


0.11 


















ELAISP90.J003244-424803 


00 


32 


44.4 


-42 


48 


3.6 


0.12 


1.2 




19.0 


0.192 










ELAISP90.J003323-432634 


00 


33 


23.2 


-43 


26 


34.8 


0.19 






18.7 


0.316 











Table 1. ELAIS SI 90^m sample. The 90^m photometry has a ~ 30% possible systematic uncertainty (see Paper III), and 
the 15/^m photometry is subject to a possible systematic scaling to fainter fluxes by up to a factor of ~ 1.5 (see Paper II). We 
assume Ho = 65 km s _1 Mpc -1 , Slo = 1, A = 0, and bolometric luminosities are calculated at 90/im assuming a constant 
vL v spectrum. The maximum accessible redshifts z max quoted assume no evolution, and are labelled respectively as P, R, 
O and C for maximum redshifts imposed by the 90^m PHOT flux limit, the radio flux limit, the optical magnitude limit 
and the CAM flux limit. Two galaxies (ELAISP90_J003721-434228 and ELAISP90.J003021-423657) are Seyfert II galaxies 
and two further galaxies (ELAISP90.J003431-433806 and ELAISP90.J003133-425100) have an early type spectrum; all 
the remainder show starburst emission line spectra. Also listed are 4 redshifts for candidate optical IDs with the optical 
magnitudes listed, but which lack other multi-wavelength detections and are too faint optically to be certain that the 
correct ID has been found. All 4 such identifications show starburst emission line spectra, but nevertheless are excluded 
from the luminosity function analysis as they fail the strict selection function defined in the text. 



than 7? = 17, or (b) a reliable 15/xm or 1.4GHz iden- 
tification, together with an optical ID brighter than 
R — 20. The 15pim identifications were made to S15 > 
2mJy, and the radio identifications to 0.2 — 0.4mJy 
depending on position (Gruppioni et al. 1999). Only 
90/im galaxies in the area with radio and 15/im survey 
coverage were considered. We use the M82 starburst 
model from Rowan-Robinson et al. 1997 to estimate 
the K-corrections, with Seo/jm = 1205*1. 4GHz and a 
radio spectral index of dlog 10 S/dlog 1Q v = —0.8. For 
optical K-corrections we assume an optical spectral 



index of —3. The luminosity-redshift plane is shown 
in figure [j], with bolometric luminosities calculated at 
90/j.m assuming constant vL v (for ease of conversion 
to other luminosity scales). Our method for calculat- 
ing the luminosity function in the face of this com- 
plicated selection function is dealt with in the ap- 
pendix. Note that our method relies on the under- 
lying assumption that no galaxies are missing at all 
redshifts due to the multivariate flux limits. However, 
we can be confident of this in our case, as any suffi- 
ciently low-redshift galaxy will have an optical iden- 
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Differential LF Differential LF (1+z) 3 ' evolution 




L(9Oyum)/L0 L(9O|xm)/L0 



Figure 2. Luminosity function assuming no evolution (left) and (l + z) 3 luminosity evolution (right). Upper figures show 
the integral luminosity function, lower figures show the differential form. Shaded area is the ±1<t error region from this 
study. For the differential counts, the shaded area indicates the error on the luminosity function in a ±0.5 dex bin; the 
differential counts are therefore effectively smoothed by a 1 dex boxcar. The data points show the local lOO^tm luminosity 
function of Rowan-Robinson, Helou & Walker (1987), and the full line is a fit to this local data. The broken lines show the 
populations in the Rowan-Robinson (2000) model. The long dashed line is the "cirrus-like" population, i.e. galaxies with 
bolometric luminosities dominated by new stars heating previously-created dust. N.B., this population is not necessarily 
identified with NGC 6090-like cirrus galaxies. The short dashed line is the starburst population; dash-dot indicates Arp- 
220-like galaxies; dash-dot-dot-dot shows the AGN population. 



tification that will pass criterion (a) of our selection 
function. A hypothetical local population of very op- 
tically faint but far-IR-bright galaxies can already be 
excluded from IRAS samples. We must also assume 
that our 89% complete optical spectroscopy is a ran- 
dom sparse sample of the total identifications (though 
not necessarily of the total 90/^m sample). Although 
with our current spectroscopic sample this is an a 
posteriori selection rather than a priori, there were 
no selection biases in the spectroscopic data acquisi- 
tion which could skew the sample selection function. 
We therefore correct our effective areal coverage by a 



factor of 0.89 to account for the small optical spec- 
troscopic incompleteness. 

The integral luminosity function is given by 

$(> l) = J2 c,, (i) 

Li>L 

where the sum is performed over all objects having 
luminosities greater than L. In figure ^ we show this 
integral luminosity function, and compare it to the 
local 100/im luminosity function derived by Rowan- 
Robinson, Helou & Walker 1987 (binned data). The 
latter was derived from a sample in the north Galactic 
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Figure 3. ELAIS source counts at 90/xm, compared to a 
no-evolution model (shallow curve) and the (1 + z) 3 lumi- 
nosity evolution model (steeper curve). The no-evolution 
curve has a slightly shallower slope than its Euclidean 
equivalent. Also shown as open symbols are the Lockman 
Hole counts from Linden- V0rnle et al. (2000). All errors 
are Poissonian, except in the case of a single object in a 
bin for which the ±1ct bounds on the number of objects 
in the bin are 0.18 — 3.3. 



pole, and is subject to possible large scale structure 
variations. To estimate this, we compared the normal- 
isation of their 60/im luminosity function with that of 
Saunders et al. (1990), and derived a correction factor 
of 0.77 to the normalisation of their 100/^m luminos- 
ity function. The best-fit function to this local data 
is plotted as a full line. The adopted functional form 
is identical to Saunders et al. (1990): 



4>{L) = faij^-) 1 a exp( 



2<r 2 



log 2 (l 



L 



)) 



(2) 



and the best-fit parameters are given in table ^. 
Also plotted in this figure are the models of Rowan- 
Robinson (2000). 

The shaded area shows the ±l<r errors from our 
sample. There is a marginally significant excess at the 
highest luminosities. The {V/Vmax) statistic for this 
sample is 0.54 ± 0.07; a Kolmogorov-Smirnov test on 
the V/Vmax distribution shows only a 34% probability 
of inconsistency with the top-hat distribution U[0, 1]. 
Assuming (1+z) 3 luminosity evolution (figure^) gives 
(VyVmax) = 0.51 ±0.07, and in both the evolving and 
non-evolving models, a Kolmogorov-Smirnov test on 
the observed luminosity distribution also gives accept- 
able confidence levels. Our spectroscopic sub-sample 
on its own is therefore not sufficiently large to reliably 
detect evolution. 

A much stronger constraint on the strength of the 
evolution comes from the source counts slope. The 
90/im counts in the entire ELAIS areas are signifi- 
cantly non-Euclidean (Paper III), but the counts do 
show a suprising upturn at the faintest end (figure 
^|). The cause of this upturn is not clear, but it is too 
large to be entirely attributable to uncertainties in the 



0.100 r 




2.0 



Figure 4. Predicted cumulative redshift distribution for 
surveys at 100k mJy (diagonal hatching), 10k mJy (hor- 
izontal hatching) and Ik mJy (vertical hatching), where 
k = 1 ± 0.3 indicates the uncertainty in the ELAIS PHOT 
absolute flux calibration. The bounds of the shaded areas 
are dominated by the assumed maximum redshift at which 
luminosity evolution occurs. 



Phi 



5.4 ±1. 



,L*hl 5 /L e 



9.67 ± 1.47 



1.73 ± 0.04 



Table 2. Best-fit parameters of the local luminosity func- 
tion, with a = 0.724 assumed. The reduced x 2 of the best 
fit is 0.66. 



completeness correction. One possibility is that some 
of the faintest sources are spurious glitch events, but 
all the sources have been eyeballed and accepted by 
at least two observers. Such a glitch population would 
have to be somewhat pathological. Evidence for com- 
parably strong evolution at < 200mJy in the Lock- 
man Hole is also reported by Matsuhara et al. (2000) 
based on model fits to the fluctuations in their 90^im 
maps, consistent with the faintest ELAIS points. Nev- 
ertheless, in the absence of a clear model-independent 
interpretation, we confine our source count fits to the 
> 200m Jy region. We note that the sources used in 
the luminosity function are all reliably cross-identified 
at other wavelengths. The source counts at > 200mJy 
are in good agreement with the Lockman Hole counts 
from Linden- V0rnle et al. 2000 (figure |^) who use 
a very different source extraction and flux calibra- 
tion procedure to ELAIS. The ELAIS and Lockman 
counts for the faintest bin in figure ^ (165 — 294mJy) 
are 27 ± 4 and 22 ± 9 deg -2 dex -1 respectively, and 
in the brightest bin (294 - 522 mJy) 19 ± 4 and 4±3. 3 
deg -2 dex -1 respectively. Note that this represents 
only a single source in the Lockman Hole, so Poisso- 
nian errors would underestimate the error in this case. 
The ±1<7 bounds in the case of 1 observed source are 
0.18 — 3.3, and all the remaining errors quoted in the 
counts are ±la Poissonian errors. 

If we assume the local luminosity function de- 
rived above undergoes (1 + z) a luminosity evolution 
we can obtain a constraint on the evolution param- 
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eter a. Note that this is relatively insensitive to the 
flux calibration uncertainty (Paper III). The correct 
relative normalisation between 90/mi and 100/im for 
these purposes is one which achieves continuity in the 
source counts. The current 90^im calibration satisfies 
this almost exactly (Paper III) in that the bright end 
of the 90^tm PHOT counts dovetail with the faint end 
of the 100/im IRAS counts, with little room for error 
in the relative calibration. The Kolmogorov-Smirnov 
test on the shape of the source counts yields a lower 
limit of a > 1.5 at 95% confidence, and no-evolution 
models can be excluded at the 98% level. Stronger 
constraints still are possible if we include the nor- 
malisation of the source counts. Our goodness-of-fit 
statistic for this is Pks x Pa, where Pks is the signif- 
icance level of rejecting the model of the shape of the 
source counts, and P^ is the significance level of re- 
jecting the predicted normalisation. The distribution 
function for such statistics is quoted in e.g. Dunlop & 
Peacock 1990. We obtain a = 2.45 ± 0.85 (95% confi- 
dence). This is stronger than the constraint obtained 
purely from a V/Vma,^ analysis of the spectroscopic 
sub-sample, partly because the number of galaxies 
comprising the source counts is much larger, partly 
also because we have assumed a fixed functional form 
for the zero redshift luminosity function, but mainly 
because we have set the z = normalisation from 
IRAS. 



4 DISCUSSION AND CONCLUSIONS 

We have found that both the shape and normalisa- 
tion of the 90/im luminosity function is well-described 
by the local luminosity function of Rowan- Robinson, 
Helou and Walker (1987), with evolution consistent 
with (1 + 2) 2 - 45 ± ' 85 p Ure luminosity evolution, using 
the transformation S100 ~ Sgo- Pure density evolution 
is already excluded in this population as it conflicts 
with sub-mjy radio source counts (Rowan-Robinson 
et al. 1993). 

We can use these results to predict the redshift 
ranges of current and future 90/im surveys. We as- 
sume the evolution extends to between z — 1 — 2 and 
is then unevolving. Figure ^| shows the 95% confidence 
limit for surveys of various depths, incorporating the 
uncertainty in the ELAIS 90^im absolute flux calibra- 
tion. One corollary is that at least 10% of the cur- 
rently unidentified 90^im galaxies in ELAIS extend to 
z > 0.5, possibly even z > 2. 

Is it reasonable to assume that star formation 
dominates the bolometric luminosities of our 90^im 
sample? Genzel et al. (1998), Lutz et al. (1998) and 
Rigopoulou et al. (1999) have presented mid-infrared 
spectroscopy of ultraluminous galaxies (L > 10 12 Lq) 
in the same redshift interval as our sample, and using 
the PAH features as a star formation indicator found 
that star formation dominates the bolometric power 
outputs in 70 - 80% of ULIRGs. The AGN bolomet- 



ric fraction derived this way decreases monotonically 
at lower luminosities (e.g. Lutz et al. 1998). In local 
IR galaxies between 10 11 and 1O 12 L0, ~ 15% have 
Seyfert II spectra (Telesco 1988), though care should 
be taken in interpreting AGN dominance from optical 
spectra (e.g. Lutz et al. 1999, Taniguchi et al. 1999). 
At L < 10 11 Lq most local IR galaxies are single, gas 
rich spirals (e.g. Sanders & Mirabel 1996) in which 
many lines of argument point to starburst dominance 
in the IR (e.g. Telesco 1988): for example, the simi- 
larity of the spectral energy distributions to those of 
galactic star forming regions; the linear LrR— Leo cor- 
relation over many orders of magnitude in luminosity; 
and optical/IR spectroscopic confirmation of star for- 
mation activity. It therefore seems likely that our sub- 
ultraluminous population should not be powered by 
active nuclei, as is also suggested by the low fraction of 
spectra in our sample (10%) with Seyfert II features. 
(This does not imply that AGN activity cannot be 
present at a weak or bolometrically-neglectable level, 
e.g. Ho, Filippenko & Sargent 1997.) At the lowest 
luminosities there may be a significant contribution 
from cirrus which is at least partly illuminated by the 
old stellar population in the galaxies (Telesco 1988, 
Morel et al. 2000), though such galaxies nevertheless 
also still obey the starburst radio-far-IR relation (e.g. 
Condon 1992). We conclude that the galaxies in our 
sample have their bolometric luminosities dominated 
by star formation. 

The Saunders et al. (1990) local luminosity den- 
sity implies a local star formation rate of 0.023L0.002 
M yr" 1 Mpc" 3 (Oliver et al. 1998) assuming a 
Salpeter initial mass function from 0.1 Mq to 125 
Mq. Our results imply this rises to 0.036 ± 0.009 
M Q yr" 1 Mpc" 3 by z = 0.2 (95% confidence), sig- 
nificantly higher than the Tresse & Maddox (1998) 
Ha-based estimate of 0.022 ± 0.007 M yr" 1 Mpc" 3 . 
(Note that the two Seyfert II spectra in our sam- 
ple are at the brightest end of the luminosity func- 
tion where the contribution to the luminosity density 
is the smallest. The luminosity density is dominated 
by the galaxies around or below the break in table 
^|) We believe this is due to two factors: firstly, the 
Ha-based estimates are not immune to extinction ef- 
fects, even if performing Balmer decrement reddening 
corrections (e.g. Serjeant, Gruppioni & Oliver 1998); 
secondly, the Tresse & Maddox (1998) survey area is 
sufficiently small to be affected by large-scale struc- 
ture (e.g. Oliver, Gruppioni & Serjeant 1998) . ELAIS 
is immune to these problems, as the cosmic variance 
over the PHOT survey area is < 20% at z < 0.3 (Pa- 
per I). These star formation rates are however con- 
sistent with extrapolations from those derived from 
mid-IR data at higher redshift by Flores et al. (1999) 
and Rowan- Robinson (1997). 

The prospects for improving on the results pre- 
sented here are excellent. A large-scale spectroscopic 
follow-up of the ELAIS northern areas is currently un- 
derway, which will have a major impact on the science 
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analysis of ELAIS. In particular, the ELAIS 15/im 
luminosity function will probe higher redshifts than 
accessible to PHOT at 90^im, from which the 15/im 
luminosity density can constrain the cosmic star for- 
mation history (e.g. Serjeant 2000). 
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APPENDIX: 1/Vmax WITH 
MULTIVARIATE FLUX LIMITS 

In a volume-limited non-evolving sample, the comov- 
ing number density of objects can be calculated triv- 
ially: 



Vrn 



N 
V 



N 

-y 

1=1 



(3) 



where V is the volume and N the number of objects. 
Suppose that the sample is incomplete in whatever 
way, and that the probability of the i th object being 
contained in the sample is pi. Provided that none of 
the pis in the underlying galaxy population are zero, 
the comoving number density in an observed sample 
of iV b s objects can be expressed as 



(4) 



This will be an unbiased estimator of the underlying 
value. 

Applying this formalism to our current data set, 
the incompleteness is due to the multivariate flux lim- 
its. We treat the (non-evolving) galaxies as sampling 
random redshifts within the volume, so that each pi 
is the probability that such a galaxy would lie above 
the flux limits. The statement that the underlying piS 
are all non-zero is equivalent to assuming that there 
are no galaxies (in the 90/im luminosity range being 
considered) that would lie outside the selection cri- 
teria at all redshifts. We can be confident that this 
is the case, because any sufficiently local galaxy will 
have an optical identification passing criterion (a) of 
our selection function. A hypothetical population of 
FIR-luminous, optically-faint galaxies can already be 
excluded from IRAS. 

There are two equivalent ways of calculating the 
PiS. If we embed our flux-limited sample in a larger 
volume of size Vo, we can use the assumption that the 
population is non-evolving and the fact that the flux 
limits are monotonic in z to express the piS as 



V 



(5) 



where V max ,i is the volume enclosed by the maximum 
redshift z max ,i at which the i th object is visible. Note 
that this is the smallest redshift at which the object 
fails any of the selection criteria. The number density 
in a sample of iVobs galaxies is then 



«obs iV ob « 



Knax.z Knax. . 

i=l i=l 



(6) 



The RMS error is simply 



A<)> : 



\ 



V 2 



(7) 



An alternative but equivalent method of calculating 
the is to model the incompleteness due to some 
or all of the multivariate flux limits by introducing a 
weighting factor to the differential volume elements: 



dV' = 7(s, ...)dV. 



(8) 



The max ,i values would then be calculated using the 
remaining un-modelled flux limits only. If all the flux 
limits have been modelled, then Vi z msx ,i = zo where 
V(zo) = Vo- The 7 factor would drop to zero at some 
redshift if (under the previous method for calculating 
the z max ,i) zo > max(z maXi i). 

Here, our approach is to use the first method 
(equation |^) to treat the multi-variate flux limits, as 
it is less model-dependent. However we use a weight- 
ing function on the volume elements to correct for the 
ELAIS 90/^m completeness function (figure 6 of Paper 
III): 



V' ■ = 

r in ax. j, 



7(; 



, bi)——dz 

dz 



(9) 



where Si is the 90/xm flux of the i* object. 

To calculate the luminosity function in a redshift 
bin, an additional top-hat selection function is applied 
in redshift space and incorporated into the calculation 
of the z max ,i. The volume enclosed by the minimum 
redshift of the bin must also be subtracted from the 
Vrnax.i in equation [| We are currently confident that 
none of the underlying piS are zero, but if we cal- 
culated the luminosity function in a series of redshift 
bins this would not necessarily be the case for all bins. 
For example, there are galaxies in table |l| which only 
pass criterion (a) at the lowest redshifts, and these 
may be unobservable in principle in higher redshift 
bins. We would then need to make some model of 
the multi-variate correlations to correct for the miss- 
ing pi=0 galaxies. We therefore restrict ourselves to a 
single redshift range (de-evolving the galaxies if nec- 
essary) , partly to avoid the model-dependence of cor- 
recting for missing populations, but partly also be- 
cause of our small sample size. 
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